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Abstract 

We study the dynamics of a spherically symmetric thin shell of perfect fluid em¬ 
bedded in d— dimensional Anti-de Sitter space-time. In global coordinates, besides 
collapsing solutions, oscillating solutions are found where the shell bounces back and 
forth between two radii. The parameter space where these oscillating solutions exist is 
scanned in arbitrary number of dimensions. As expected AdS 3 appears to be singled 
out. 


1 Introduction 

The non-linearity inherent to Einstein field equations makes the obtention of exact analytical 
solutions a hazardous program. Among the simplifications that allow to ease this task, the 
construction of thin shell space-times stands out as a manageable approximation. Although 
idealised, thin shells are extremely useful, both from a conceptual and a computational 
perspective, since they provide nice solvable models in which trademark processes, such as 
gravitational collapse, can be explicitly addressed [I]. 

On the other hand, solving gravitational problems in Anti-de Sitter space (AdS) has be¬ 
come a theoretical laboratory were one may want to test results and solutions known from 
asymptotically flat space-times. Needless to say, this interested has been triggered by the 
so called AdS/CFT correspondence, also known as holographic gauge/gravity duality |2|. In 
such a context, a static black hole in AdS is dual to a thermal state in the CFT that lives in 
the AdS boundary. Pushing this identification beyond the static situation, it is natural to 
conjecture that thermal physics out of equilibrium in CFT can be modelled by dynamical 
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gravitational processes that involve time dependent horizons. For small perturbations this 
has led to successful computation of hydrodynamic transport coefficients from quasi-normal 
modes |3|. Large gravitational collapses of some matter field configuration should provide 
a dual description of the thermalisation process of a strongly coupled CFT from some out- 
of-equilibrium initial state [4|. This relaxation process is inaccessible to perturbative field 
theory techniques, but it is also difficult to treat from the gravitational point of view due to 
the non-linearity of the partial differential equations involved. 


This is the reason why collapsing thin shells in asymptotically AdS spaces have been widely 
studied in the framework of gauge/gravity duality. The simplest examples involve flat shells 
that fall inside planar AdSrf +2 space, where the dual quantum field theory lives on the (d+1)- 
dimensional Minkowskian boundary. Concerning the kind of matter that makes up the shell, 
null dust leading to Vaidya space-times has been considered (5, 6|. More general kinds of 
matter have also been treated recently |7|. 


Collapsing shells in global AdS should model holographically the relaxation of an isolated 
quantum system of finite size. Having a boundary conformal to Sj x M, the natural ex¬ 
pectation is that new phenomena may arise, as the sphere introduces a new length scale 
into the game. This expectation is reinforced by the fact that, with the help of numerical 
techniques, exotic behaviours have been observed in collapses involving a massless scalar 
field. In the Poincare patch, both a flat shell and a massless scalar pulse lead to direct black 
hole formation after the first infall towards the Poincare horizon |8|. In global coordinates, 
however, the same matter configuration may exhibit a transient oscillatory behaviour if its 
mass is sufficiently small. The thick shell of massless scalar field bounces back and forth 
between the origin and the AdS boundary several times while it thins out, before ending up 
forming a black hole [9,10. The initial simulations led the authors of Ref. |9] to conjecture 
that AdS is non-linearly unstable towards formation of a black hole no matter how small the 
initial perturbation is. Later on, some numerical evidence was provided to believe that there 
are “islands of stability” in the space of initial conditions. The situation at present is not 
settled but there are evidences that suggest that these islands contain each one an exactly 
periodic solution which governs the stability. In some cases, it has been posiblc to construct 
explicitly 11 such exactly periodic solutions. 


Our initial intention was to look for a simple analytical model that could reproduce such 
periodic behaviour. This led us to consider spherically symmetric d —dimensional thin shells 
embedded in global AdS ^+2 space-time. Restricting to a family of linear equations of state, 
we will show that there are regions in parameter space where the shell undergoes an exactly 
periodic motion. It is worth stressing that, in contrast, in planar AdS shells never bounce 
back once infalling 
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2 Shell Dynamics 


The shell world-volume £ is a codimension-1 hypersurface that divides the (d+2) - di mens ion al 
background space-time M. in two distinct regions: outside, _M + , and inside, Ad_. Due to 
the spherical symmetry of the problem, we know, by Birkhoff’s theorem, that the space¬ 
time metric g takes the Schwarzschild-AdS form on both Ad+, Ad_. Choosing standard 
Schwarzschild coordinates x± = (t±,r±,9i,...,9a ) to cover M.±, we find that, in this partic¬ 
ular coordinate system 


ds 2 ± = —f±(r±)dt± + f±(r±) 1 dr 2 ± + r 2 ± cin 2 d ( 1 ) 


where 


f ±( r ±) 


1 + A 


m± 

2 

r ± 


( 2 ) 


As usual, dfl d is the metric of a unit round d-dimensional sphere, and the AdS radius l is 
related to the cosmological constant A = — . In what follows, we will restrict ourselves 

to the case where the shell has no influence on the cosmological constant, so we are going 
to set l + = Z_ = l = 1 by an appropriate choice of units. Furthermore, we assume that the 
space-time inside the shell is empty AdSd + 2 and hence fix m_ = 0. In such case, m + = m 
sets the total ADM mass of the system. 


Let the shell world-volume £ be parameterized with coordinates y = (r,0i,...,0d), where 
r is the proper time of an comoving observer. The shell embedding in the ambient space- 
time is given parametrically by the function 

x s (y) = (Cfc.sOr), t ±, s (t), 9 u ..., 9 d ) . (3) 

The tangent space T P A4 of any point p G £ admits a basis formed by d+1 vectors e a = e“<9 x a, 
tangent to £, and one vector n = n a d x a , orthogonal to £. Explicitly, 

6t,± i±,s&t± T (4) 

e 6i,± = 9 g . (5) 

n± = + {f±y±, a d t ± + f±, a i±, a d r± ) (6) 

The overall positive sign of n± is fixed by requiring that n± is always directed from Ad- to 

M + Q. 

The embedding ([3]) is not arbitrary: in order for the whole space-time Ai to solve Ein¬ 
stein equations, the so called Israel junction conditions must be satisfied (see Ref. [l]). The 
first junction condition states that the induced metric h ag on £ must be continuous across £ 

[hab] = 0 (7) 
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where the brackets stand for the jump. The second junction condition relates the jump of 
the extrinsic curvature K ab with the matter composition of the shell, 


[K ab - h ab K] = -8irGS ab = -S ab ( 8 ) 

where K = h ab K ab , S ab is the shell energy-momentum tensor and we have chosen units such 
that 87 tG = 1. Projecting g onto S to find the induced metric h ab = g a /3 e< g e b we get 

dh± = h± ab dy a dy b = (-f±J 2 ±, s + f±yL,) dr 2 + r 2 ±s dQ 2 d . (9) 


The choice of r as comoving time fixes h TT = — 1, whence it follows that 

At 




f±,s 


with 


At - yJ7±,. 


+ rl s 


( 10 ) 


( 11 ) 


We have taken the positive root of t± jS , as we want the shell trajectory to be future oriented. 
Equation (10) together ( fTT| ) accomplishes two tasks. It ensures that the rr component of 
the first junction condition (|7j) is satisfied, and gives the correct normalisation to the vector 
n in 0, n 2 = 1. It also implies that it is imposiblc to cover the entire space-time M. with a 
globally defined time-like Schwarzschild coordinate, as the embedding functions t± tS (r) will 
differ at the shell. On the other hand, the radial coordinate r± has to be continuous since 
r +]S (r) = r_ ]S (r) = r s (r) must hold to signal unambiguously the shell’s radial position. This 
condition, together with ensures that all components of ([ 7 ]) are satisfied. From now 
on, we take these facts into account and change correspondingly our Ai± coordinate system 
to x± = (t±,r,9i). 


The extrisic curvature is the pullback of the Lie derivative of the ambient metric g along n. 
Several equivalent expressions can be found in the literature [I] 


Kab ~ 2 e « e 6 (^ng) a ,p ~ 

where the orthogonality condition e“n 
components and trace are 


e a e fV cv n / 3 = -rif. 


dyadyb + af) a b 


= 0 is used. In our particular setup ([9]), its 


( 12 ) 

non-zero 


K T - — 
K± ’ t ~ G 


K di - — 

1 S 


K = — + d— 
r s r s 


(13) 


The diagonal nature of Kg, together with the second Israel junction condition (J8|) , leave little 
room for the form of the shell stress-energy tensor Sg, which must be of the perfect fluid 
form 


Sg = diag 


(14) 


4 






where a will be the shell energy density and p the shell pressure. Due to spherical symmetry, 
p is independent of the particular angular direction considered. In components, ([8]) reads 
now 


[Ki\ 

IK] 




d — 1 

p + KT a 


(15) 

(16) 


With our choices for f±, we always have /3 + < (3- and, therefore, a > 0. As usual, equations 
(15), (16) need to be supplemented with an equation of state which relates the shell energy 
density and pressure. At this point, we introduce a simplification by restricting our analysis 
to the case where this equation of state is linear. In AdSrf +2 we shall write 


a 

p= d a ‘ 


(17) 


The parameter a determines the kind of matter the shell is made of. Taking a G [0, 1] it 
interpolates between dust (a = 0) and conformal matter (a = 1). This choice of equation 
of state, together with the positivity of a, implies that a + p > 0, so that the weak energy 
condition is respected. With the choice ( [l7| , equations (15), (16) can be solved explicitly. 
The final result is that the shell dynamics is fully equivalent to the one-dimensional motion 
of a particle in an effective potential 14 // 


r, + V./I = 0 


(18) 


where 

14 // = 1 + r 2 s - \mr l ~ d - - \ M 2 r :^ d ~ 1+a) . 

e// s 2 AM 2 4 

M is an integration constant that sets the shell’s proper energy E, defined as E = 
vol{S ( j)dr~ a M , where vol(Sd) is the volume of the d-dimensional unit sphere. 


(19) 

vol(S d )r d ,o = 


Notice that the potential 14 // is invariant under m —> m, M —> || and a —> — (d — 1 + a) 
so any result we may obtain is also going to hold in the range a G [— d, —d + 1], modulo the 
appropriate M redefinition. The weak energy condition will be still satisfied. 


3 Oscillating Solutions 

From ( [18] ), we know that the region where the shell is allowed to move is the -possibly 
disconnected- set of radial intervals for which 14 // < 0. The asymptotic behaviour of 14 // 
is as follows: 


• K// —> —oo as r — > 0, with 14// ~ — l/4M 2 r 2 ^ d 1+Q i for r <C 1 
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• K// —* oo as r —* oo, with Vj,// ~ r 2 for r 1 

Oscillating shell trajectories, if they exist, are confined to an intermediate radial region 
where the potential V e ff develops a well. The turning points are given by two radii r± where 
K//(' r ±) — 0 (subindices here do not refer to the inner or outer regions to the shell). The 
goal now is, fixing d, a and the space-time ADM mass m, find in what range of M oscillating 
solutions appear. It turns out that this M region is bounded by two shell rest energies M i ui 
with Mi < M u , such that 

• For M = M u there is a local minimum of V e ff touching the V e ff = 0 axis. This 
corresponds to a shell in equilibrium. 

• For M = Mi there is a local maximum of V e ff touching the V e ff = 0 axis. This 
corresponds to the transition between oscillatory and collapsing behaviour. 


Veff 



Figure 1: Typical form of the potential V e ff in the limiting cases that bound the oscillating 
shell existence region. In black, the potential at M = M u . In gray, the potential at M = Mi. 
Dashed gray lines correspond to ri and r u . The parameter values are d = 3, a = 0.5, m = 0.1, 
M u = 0.0359, Mi = 0.0346, n = 0.426, r u = 0.901. 


Fig. [l] depicts both limiting situations on AdS, 5 . To find the values of M iu , we have to 
solve the system of equations given by 14// = d r V e fj = 0. The solution is easily obtained in 
implicit form by taking m = m(d,a,r), M = M(d,a,r), which will be called the existence 
curves 


m(d, a , r) 


M(d, a, r) 


Ar d 1 (q; — (1 — a)r 2 )(d — 1 + a + (d + a)r 2 ) 
(d - 1 + 2a) 2 (1 + r 2 ) 

2r d ~ 1+a (a — (1 — a)r 2 ) 

(d — 1 + 2a)\/l + r 2 


( 20 ) 

( 21 ) 


As the d = 1 existence curves display peculiar properties we will discuss this case separately. 
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3.1 Oscillating shells in d > 1 


The existence curves allow a straightforward computation of Mi and M u . The procedure 
is illustrated in Fig. [2] (left). First, we choose some m and solve numerically the equation 
m = m(d, a, r ) for r at fixed d and a. The output are two radii, and r u , such that r; < r u : 
r u signals the position of the axis-touching minimum of V e ff(d,a,m,M u ), while r'i signals 
the position of the axis-touching maximum of V e ff(d, a,m, Mi) -see Fig. [I]-. Inserting r/ >u 
into equation (21) gives back the numerical values of Mi >u and fixes completely the form the 
potential V e ff. For any M G (Mi,M u ), it is guaranteed that V e ff possesses an oscillating 
solution. 


Looking at Fig. [2] (left), it is neatly seen that, at fixed d and a, there is a maximum mass, 
( d,a ), above which the construction just described can not be performed. Therefore, 


m. 


oscillating shell trajectories only exist for sufficiently light shells; above m max (d, a) there are 
only collapsing solutions. The fact that there are no oscillating geometries when the mass of 
the system surpasses a certain threshold is a feature that this simple model shares with more 
realistic setups, for instance, a minimally coupled massless scalar field in global AdS [9,10 


m, M M 




Figure 2: Left: existence curves for oscillating solutions in d — 3 and a = 0.99. In solid 
black, m(d,a,r). In dashed black, M(d,a,r). The gray dashed curves illustrate the construc¬ 
tion described in the text. Notice that the range of M is so tiny that the two lines depicting 
Mi and M u look coincident. Right: M(3,a,m) fora = 0.2, 0.5, 0.9 -gray, black dashed and 
black curves respectively-. Any point within the region bounded by each curve corresponds to 
a particular oscillating solution. All quantities are measured in units of m max (d, a). 

Concerning M, it is important to stress that its value sets an energy scale that is independent 
of m. Tuning M deforms V e fj and, in particular, shifts the upper turning point r + . Hence, 
one can think of this parameter in terms of the initial radius where the shell is released from 
rest and starts falling and, in consequence, M would be related to its potential energy. 
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In Fig. [2] (right) we offer several parametric plots of M(d, ci, r) as function of m(d,a,r). 
In each curve, the upper branch corresponds to M u , while the lower branch corresponds to 
Mi. I 11 the m —> 0 limit, the upper M u branch asymptotes to a line, 

M u (a,m ) ~ ^ (1 — a)~ a a ^ 2 m . (22) 

Oscillating solutions exist for M in a narrow window around this upper branch of static 
shells. This window closes at m max (d,a). If the two scales set by m and M were not inde¬ 
pendent, this finite-size region where oscillating solutions reside would degenerate into a line 
and they would cease to exist. 


To start scanning the behavior of the existence region with respect to a, let us note that 
both m(d, a, r ) and M(d, a, r ) attain their maxima at the same radius 


T max {d, Cl) 


X \J 2 a(d — 1 + a) — (d + 1) + y/(d + l) 2 — 4(d — l)a — 4a 2 ) 


V2 


y/(l — Cl) (d + Cl) 


(23) 


which also controls m max (d,a), 

TTlm.ax(,di a) 4r max(d-i ®) 


d+l \/(d+ l) 2 — 4 (d — l)a — 4a 2 — (d — 1) 


(d—l + 2a) 2 


(24) 


In the conformal fluid limit a —» 1, r max (d,a ) diverges as r max (d,a) ~ f=y so we 

expect m max (d,a ) to be also divergent. In fact 


^ ^ ^rnax [ d. CX / 


l+d 

1 (!-„) 


d 2 ~ 1 \ d + 1 


1 — d 
2 


(25) 


This result shows that the maximum mass for which oscillating solutions exist grows un¬ 
bounded as ci —» 1; for d > 1, there are oscillating solutions with arbitrary high ADM 
energy as long as the shell matter is sufficiently near conformallity. In the opposite limit 
(pressure-less dust) a —» 0 


TWmax^d, Cl) ~ 


8 


d — 1 
d + l 


l+d 

2 


l+d 

a 2 


(26) 


and the allowed region for oscillating trajectories shrinks down to zero. Physically this 
means that there are no pressure-less oscillating solutions, i.e. in order to be stable against 
gravitational collapse, the shell matter must have some self-interaction. As an aside, note 
that the limiting behaviour of M(d, a, r) as a —» 0,1 just follows the behaviour of m(d, a, r ), 
since 

a (d — 1 + 2a)\f\Mr 2 
r) “ r 2(d — 1 + ci + (d + a)r 2 ) 



m(d, a, r ) . 


(27) 
















Ail important consistency check is to verify that lies outside the position of the event 
horizon, rv Otherwise, instead of having an oscillating solution, the shell trajectory would 
represent direct gravitational collapse starting from r + . As r_ > r;, it is sufficient to show 
that ri > r'h- The static event horizon location, rh(d,m), is the solution of the equation 
1 + — rnr\~ d = 0. Since we don’t know explicitly r; = ri(d,a,m), but instead m = 

m(d, a , ri), we are going to define correspondingly rrih(d, r ) = r d_1 (l + r 2 ). It is easy to see 
that, if m h fd, r ) —m(d, a, r) > 0 for all d, a, r, the consistency condition t'i > r h always holds 
-see Fig. [3}. We find that 


m h (d, r ) 


m(d, a, r ) 


r d -i (d — l + (d+ 1 )r 2 Y 
(1 + r 2 ) (d - 1 + 2a) 2 


(28) 


which is positive definite. Thus, every oscillating shell trajectory found lies entirely outside 
the event horizon of the would-be black hole. 


m, m.h 



Figure 3: Position of the event horizon -gray- versus m existence curve -black- in d = 3 for 
a = 0.9. The line m = 1.5 is drawn -gray dashed-, as well as its intersection points with m, 
ri and rrih, rh -gray dashed- 


3.2 Oscillating shells in d = 1 

The case d — 1 is special because, as mentioned before, some general properties seen in d > 1 
do not hold anymore. Representative existence curves are plotted in Fig. [4} The analytic 
form of m(d, a, r),M(d, a, r) for d — 1 is 

^4 

m(l,a,r) = m(a,r) = 1 + r 2 - ——(29) 

(1 + r z )a z 

M(l,a,r) = M(a,r) = ii(30) 

ayl + r A 
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m, M 


M 




Figure 4: Left: existence curves for oscillating solutions in d — 1 and a = 0.99. In solid 
black, m as a function of r. In dashed black, M. Right: M(l,a,m) fora = 0.2, 0.5, 0.9 
-gray, black dashed and black curves respectively-. 


It is easy to prove that ri > r^, as equation (28) still holds. The major difference with 


respect to the d > 1 case comes when we evaluating the allowed mass range for oscillating 
shell trajectories. At d — 1 


fa) = 


1 + y/1 - a 2 


cr 


(31) 


is still divergent as a —>■ 1, although in a milder way, r max (a ) ~ (1 — a) * instead of the 

(32) 


(1 — a) 2 divergence in (23). However, 

^Imaxfpf 


l + vT^ 


cr 


so, unlike the d > 1 case, for d = 1 the maximum allowed mass of an oscillating solution 
does not grow without bound as a —* 1. Instead, it goes to a finite value, m = 2 (in our 
conventions). Note also that, in the a —> 0 limit, the m-range does not close down. Instead, 
m mai (a) —» 1 as a —)■ 0. In AdS 3 there are oscillating shell solutions for any m < 1, even for 
dust made shells. 


Looking at figure Fig. [4| we observe that there is a qualitative difference between shells 
of mass above and below the threshold m — 1. Each mass m belonging to the interval 
[1, m max (a)] has two associated radii, ri and r u , such that, like in higher dimensions, ri sig¬ 
nals the position of the 14// = 0 axis-touching maximum, while r u signals the position of the 
axis-touching minimum. However, for m < 1 the axis-touching maximum disappears. This 
does not mean that shells with m < 1 are allowed to reach the point r = 0, because the po¬ 
tential barrier does not vanish at any finite M: as M —)■ 0, the maximum of the barrier gets 
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radially displaced towards the origin and tends to a constant value. At the same time, r + 
grows unbounded. In the same way that oscillating shells can not cross their Schwarzschild 
radius, in AdSs there is also a mechanism forbidding the possibility of reaching r = 0 : the 
shells can not form naked singularities. Again this is in parallel with the phenomenology 
shown by a massless scalar pulse in AdSs [l3,14|. The behaviour of the barrier is depicted 
in figure Fig. [5] 


It is easy to explain the behaviour of V e ff as M —> 0 in analytic terms. Let rbarrier be the 
position of the maximum of the potential barrier. If r « 1, d r V eff ~ a (M 2 r- 2 “ - ^r 2 “), 


so 


r ba 


M 


Therefore, Terrier —* 0 as M —> 0. At r = r barrier , the potential is finite 

F barrier = f{r barrier) ~ 1 rn . 


(33) 


(34) 


Regarding the large r turning point, r + , note that if V has a root r + such that r + ^ 1, we 
have that V (r + ) ~ r 2 — jrAr 2 " = 0 , which is solved by 


r 4 - = 


m \ i -a 


2 M 


(35) 


This proves that r + —> oo as M —> 0. 


VeH Vsff 



Figure 5: Left: for a = 0.99 and m = 0.8, potential V e ff at M = 0.3899, which corresponds 
to a static shell located at r u = 7.699. Right: evolution of the potential barrier for the 
oscillating shell at different M. The black curve corresponds to V e ff on the left picture. 
From left to right, gray curves represent potentials with M = 0.3, M = 0.2, M = 0.1 and 
M = 0.05. Note that r bar rier gets smaller as M does, while Vb arr ier —» 1 — m = 0.2 in this 
case -black dashed- 
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4 Conclusions 


We have explored a family of solutions of General Relativity in global AdS exhibiting periodic 
behaviour. It involves thin spherical shells of matter governed by linear equations of state 
that range from dust to conformal fluids. We have provided exact analytic expressions for 
the regions of existence of such solutions and analysed their behaviour in certain limits. The 
conclusions match in a natural way with similar results obtained in numerical simulations of 
realistic shells of massless scalar held. In particular, we have checked that oscillating shells 
never reach the Schwarzschild horizon. The case of AdS 3 has been separately analysed due 
to the presence of the mass gap. 

As several simplifying assumptions have been made, each one can be relaxed independently 
in order to test the robustness of the results obtained here. For instance, more general 
equations of state, i.e. polytropes, can be considered. The requirement that the interior 
space-time J\4_ is an empty AdS space with the same cosmological constant as the exterior 
one can also be lifted. Finally going beyond, spherical symmetry is a natural extension of 
this work. While the case for rotating oscillating solutions looks feasible, the addition of 
angular deformations appears as a formidable challenge. 

In the light of gauge/gravity duality, these oscillating geometries should correspond to out- 
of-equilibrium field theory states that never thermalise [l5]. As has been stressed, the simple 
nature of the solutions found would allow a relatively easy calculation of holographic proxies 
of field theory quantities (such as entanglement entropy |7j, two-point correlation functions 
or mutual information) that would help to characterise the non-thermalising state. This 
issue is currently under research. 
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